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ABSTRACT 

We present a photometric investigation of the globular cluster population of NGC 
1399, the central galaxy in the Fornax cluster, in Washington C and Kron R filters within 
a field of 36' x 36', corresponding to about 200x200 kpc at the Fornax distance. This 
is the largest area around this galaxy ever studied with CCD photometry. The cluster 
system of NGC 1399 is found to extend further than 100 kpc away from the galaxy. The 
color distribution exhibits a pronounced bimodality. Within a radial distance of about 
55 kpc, the red clusters are more strongly concentrated to the center than the blue 
clusters. At larger radii, the surface density profiles of the clusters are indistinguishable 
and match well the distribution of the galaxy light. Over the entire radial distance 
range, the surface brightness profile of NGC 1399 can be very well fitted by a power law 
with an exponent of -1.85 and a core radius of 3.3". No steepening of the luminosity 
profile can be detected at large radii. We suggest that the power law profile of NGC 1399 
results from the galaxy being embedded in a large dark matter halo, which prevents the 
stellar density profile from steepening outwards. The cluster system contains 6450 it 700 
clusters and the specific frequency is found to be 5.1 it 1.2 in the V band. While NGC 
1399 shows a pronounced color gradient the nearby comparison galaxy NGC 1404 does 
not show such a gradient. Using simple assumptions about the underlying population 
that formed during the same star formation event as the globular clusters, we present 
a model in which we use radially changing local specific frequencies for the red and 
blue subpopulations to fit the observations. We find that within 7' the required specific 
frequency of the blue clusters alone is a factor of approximately 3 larger than that of 
the red ones. Outside this radius however, both populations have a same high local 
specific frequency of around 8 and 13 (blue and red clusters, respectively). 

Subject headings: galaxies: individual (NGC 1399, NGC 1404) — galaxies: elliptical 
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and lenticular, cD — galaxies: clusters: general — galaxies: photometry — galaxies: 
halos 

1. Introduction 

One of the central and challenging topics in the context of the structure and evolution of 
early-type galaxies is the investigation of their globular clusters systems (GCSs), how they formed, 
and what they might tell us about the formation of their host galaxies (for recent reviews on the 
subject, see Ashman & Zepf (1998), Elmegreen (1999), van den Bergh (2000), Harris (2001)). 

It has now emerged from numerous studies (e.g. Kundu k. Whitmore (2001a), Kundu & 
Whitmore (2001b), Larsen et al.(2001), Gebhardt & Kissler-Patig (1999), Geisler et al. (1996), 
Whitmore et al. (1995), Zepf &; Ashman (1993)) that GCSs of early-type galaxies are by no means 
homogeneous populations, but exhibit substructure in terms of metallicity, kinematics and spatial 
distribution. The formation of massive clusters does not occur in isolation, but probably even needs 
to be accompanied by strong star forming activity in general. This is apparent in ongoing merger 
events (see the list in Harris (2001)) but also in the case of "normal" spiral galaxies (Larsen &: 
Richtler (1999), (2000)). Therefore, the presence of GCs must be indicative of a corresponding field 
population with similar ages and metallicities and thus GCs might be considered as tracers of this 
field population, which otherwise is unaccessible due to the low surface brightness of galaxy halos 
and/or because it is mixed up with other populations. 

It might then be possible to trace back the population composition of an elliptical galaxy, the 
distinct metallicities and spatial distributions of the subcomponents, and compare these properties 
to what is expected if, for example, the elliptical galaxy is a result of a merger (or multiple merger) 
of disk galaxies (Bekki et al. (2002)). However, the relation of a GCS and its components to the 
field components can be complicated. Even if the efficiency of cluster formation per unit mass had 
a universal value, as argued by McLaughlin ((1999)), it does not mean that one can infer directly 
the luminosity of the underlying field population from the number of globular clusters because 
generally the fate of the gas is unknown. 

A promising target to gain better insight into the connection between the relation of cluster 
subpopulations to the underlying field population is NGC 1399, the central galaxy of the Fornax 
cluster. It has long been known to possess a large number of globular clusters (Dawc &: Dickens 
(1976)) out to large radii (Hanes &: Harris (1986), Harris & Hanes (1987)). Photographic studies of 
NGC 1399 showed that its luminosity profile cannot be described by a de Vaucouleurs profile but 
by a power-law (Schombert (1986), Killeen & Bicknell (1988), Caon et al. (1994)). Whether NGC 
1399 is a cD galaxy or not depends on the applied definition, which in turn differs from author to 
author. For example. Bridges et al. (1991) classify it as a cD, while Schombert (1986) does not 
include it in his list of cD galaxies. The spatial distribution of galaxy light and globular clusters 
have been compared by Bridges et al. ((1991)) and Wagner et al. ((1991)) who performed the first 



-3- 



CCD studies of the GCS of NGC 1399. Except for the very inner regions, both studies found the 
surface density profile of the clusters to agree well with the galaxy light. Wagner et al. and Bridges 
et al. also obtained the first color information of the GCS using the B and V bands. While Wagner 
et al. could not see a color gradient, Bridges et al. detected a shift of 0.2 mag to the blue between 
0.5' and 3'. The latter authors moreover derived a turn-over magnitude of V=23.8 and a distance 
of 18.5 Mpc, which still is in very good agreement with modern values (Richtler et al. (2000)). The 
specific frequency (which is the total number of globular clusters normalized to the host galaxy's 
luminosity, Harris & Van den Bergh ((1981))), derived by Harris &Hanes ((1987)), Bridges et al. 
((1991)) and Wagner et al. ((1991)) was unusually large for an elliptical galaxy (16 it 4), a property 
which is shared with other central galaxies in galaxy clusters (e.g. Harris (2001)). 

The next contribution came from the Washington system: Ostrov et al. ((1993)) measured 
Washington colors of NGC 1399 globular clusters and also detected the color gradient reported 
earlier. They also mentioned a "multimodality" in the metallicity distribution, indicating that 
the GCS could be composed of different populations with distinct properties. Kissler-Patig et al. 
((1997)) also noted the broad V-I color distribution of the NGC 1399 clusters in comparison with 
other Fornax ellipticals and a high specific frequency, but could not confirm the color gradient seen 
by Ostrov et al. ((1993)) and Bridges et al. ((1991)). 

An improved photometric analysis in the Washington system (Ostrov et al. (1998)) showed 
that the color distribution of the clusters was rather "bi- modal" than "multi-modal". In this work, 
a rather low specific frequency of only 5.6 was derived. 

The open questions are still: How far docs the NGC 1399 GCS extend? What substructure 
can be identified and how does it relate to the host galaxy? Does the population structure change 
in the outer parts? How does the cluster distribution compare to the stellar distribution? What 
significance does the quoting of a general specific frequency have, if the GCS is substructured? 

From the photometric side, progress can only be expected if a much larger field is investigated 
than has been done in the past. With the advent of the new MOSAIC cameras, we now have the 
opportunity to do this. Moreover, the excellent metallicity resolution of the Washington photo- 
metric system enables us to reconsider the metallicity distribution based on much better number 
statistics than in previous studies. 

This paper is part of a larger effort to study NGC 1399 and its cluster population. The second 
paper in this series by Richtler et al. ((2002)) (Paper ll) will use spectroscopic observations of 
globular clusters for a dynamical study of NGC 1399. 
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2. Observation &z Reduction 

2.1. The Data and Basic Reduction 

Data from 5 observing runs are contributing to the present photometric study. The main data 
set consists of Washington wide-field images of NGC 1399 (plus a background field) taken with 
the MOSAIC camera mounted at the prime focus of the CTIO 4m Blanco telescope. Additionally, 
three sub-fields were observed at the VLT (F0RS2) in V and I. The location of the MOSAIC field 
around NGC 1399 and the VLT fields are shown in Fig. 1 as overlays on a DSS image. 

We used the Kron-Cousins R and Washington C filters, although the genuine Washington 
system uses Tl instead of R. However, Geisler (1996) has shown that the Kron-Cousins R filter 
is more efficient than Tl, due to its larger bandwidth and higher throughput, and that R and 
Tl magnitudes are closely related, with only a very small color term and zero-point difference 
(R — Tl ~ —0.02). The MOSAIC wide-field camera images a field of 36' x 36', at a pixel scale of 
0.27"/pixel. The first three observing runs were performed in the 8 channel read-out mode; for the 
last one (November 2001) the 16 channel mode was available. Further information on the MOSAIC 
camera can be found in the MOSAIC homepage ^. 

During the first run (December 1999), we observed one MOSAIC field centered near NGC 1399 
(see Fig. 1) in the C filter, while the R band images were kindly obtained by D.Zurek and M.Shara. 

In the second run (December 2000), we observed a background field 3.5° north-east of NGC 
1399 (RA= 3''48' Dec= -32°32') in C and R, as well as Washington standard fields. 

In the third run (November 2001), we observed adjacent fields to NGC 1399 that overlap with 
the central field. We also observed standard fields under photometric conditions. This run was 
used for the photometric calibration of the presented Washington data. 

In addition to the MOSAIC data, we analyzed three fields (6.8' x 6.8') near NGC 1399 for 
which 5 min exposures in V and I were taken with VLT /F0RS2 in December 2000, in the course 
of a campaign aimed at measuring radial velocities of globular clusters around NGC 1399 (Richtler 
et al. (2001)). 

In Tab. 10 we summarize the relevant information on the different observing runs. 

The MOSAIC data were reduced using the mscred package within IRAF. In particular this 
software is able to correct for the variable pixel scale across the CCD which would cause otherwise a 
4% variability of the brightness of stellar-like objects from the center to the corners. The flatfielding 
resulted in images that had remaining sensitivity variations < 1.5%. In particular Chip ^4 and 
Chip #5 showed discernible remaining flatfield structure (but within the given deviation). 

To facilitate the search for point sources on the MOSAIC images model light profiles of 



''http: / /www.noao.edu/kpno/mosaic/mosaic.html 
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NGC 1399 as well as of the neighboring galaxies NGC 1404 and NGC 1387 were subtracted. This 
was done using a median filter with an inner radius of 9.5" and an outer radius of 11". This size 
is large enough for not altering the point source photometry which has been verified with artificial 
star tests described later on. 



2.2. Identification of Point Sources 

We used SExtractor (Bertin & Arnouts (1996)) and DAOPHOT 11 for image classification 
and photometry, respectively. Our extensive tests showed that each program was optimized for 
these respective tasks. To distinguish between point sources and extended objects, we compared 
the SExtractor stellarity indices for objects foimd both on the MOSAIC and the VLT images. 
With the aid of the superior seeing on the VLT images, we determined 0.4 as a lower limit for the 
SExtractor classifier on the MOSAIC images. 

Using this classification, the comparison with the VLT fields revealed that (90 ±5)% of the 
point sources found in the MOSAIC fields were correctly classified as point sources. Only 1.2 it 0.1 
galaxies arcmin"^ have been erroneously classified as point sources in the MOSAIC images. With 
less a strict selection criterion one could lower the number of lost point sources. However, this would 
come at the expense of additional contamination by background galaxies. Since we are interested 
in the GCS at large radii, i.e. at low surface densities, we rather prefer a lower completeness over 
an increased contamination by galaxies. 

We retained 10453 point sources while 13600 resolved objects were discarded. The brightest 
non-saturated objects have Rsa 18 (slightly depending on the individual MOSAIC chip). 

For further study we chose a faint magnitude cut-off, depending on color. The limit also 
accounts for the fact that the background field is not as deep as the NGC 1399 field. We therefore 
defined this cut-off by the completeness limit of the background field and regarded only objects 
brighter than this limit (see Fig. 3) . 

2.3. Photometric Calibration 

Standard fields for the photometric calibration have been observed during the November 2001 
run. In each of the 3 nights, we observed 4-5 fields, each containing about 10 standard stars from 
the list of Geisler ((1996)), with a large coverage of airmasses (typically from 1.0 to 1.9). It was 
possible to use a single transformation for all three nights, since the coefficients derived for the 
different nights were indistinguishable within the errors. 

We derived the following relations between instrumental and standard magnitudes: 



Tl = 



Rinst + (0.72 ± 0.01) - (0.08 ± 0.01)Xr 
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+(0.021 ± 0.004) (C-Tl) 
C = Cinst - (0.06 ± 0.02) - (0.30 ± 0.01)Xc 
+(0.074 ± 0.004) (C -Tl) 

The standard deviation of the difference between our caUbrated and tabulated magnitudes is 
0.018 mag in Tl and 0.027 mag in C (see Fig. 2). 

To caHbrate the NGC 1399 field we identified isolated stars which have been observed on three 
overlapping fields during the November 2001 run. These stars were used to determine the zero 
points, while the color terms were adopted from the November 2001 run. The scatter between 
the zero points determined from individual stars is 0.03 mag and most probably due to flat field 
uncertainties. 

The final uncertainties of the zero points are 0.03 mag and 0.04 mag for R and C, respectively. 
This results in a calibration uncertainty in C-Tl of 0.05 mag (the uncertainty in the color term can 
be neglected). 

Ostrov et al. ((1998)) also obtained Washington photometry for a smaller field centered on 
NGC 1399. We used their Fig. 6 and Fig. 10 to obtain the color and luminosity distribution of 
their objects. In Fig. 10 we compare the color distribution of the same subficld, and in Fig. 9 
the luminosity functions are compared. These figures show the good agreement between the two 
photometric data sets. 

The VLT data in V and I have been calibrated using the zero points, airmass coefficients and 
color terms published on the ESO Web page for the corresponding night ^ . We used the same zero 
point for all fields despite the fact that the gain is different for all four amplifiers, since we divided 
by a normalized flatfield and thus we had to average the gain as well (in the high gain mode) . The 
ESO calibration for the I filter includes a color term Iq oc linst — 0.07(V — R). Because we had no 
VLT observations in the R band, we used a subset of clusters for which we had VLT (V,I) and 
(MOSAIC) R observations and found the the following relation between (V-I) and (V-R) colors: 
(V-R)= -(0.10 + 0.04) + (0.47 + 0. 04)(V-I), with a scatter of cry-R = 0.04 mag. This relation is used 
for the photometric calibration where no R observations were available. This procedure introduces 
an error of 0.04 mag in the final magnitudes as one can judge from the scatter of the relation. For 
Milky Way GCs the following can be obtained with the data assembled in the Mc Master catalogue 
(Harris (1996)): (V-R)= -(0.01 + 0.01) + (0.49 + 0.01)(V-I). Thus while the slope of the relation 
is the same, the V-R color is either 0.09 mag bluer in NGC 1399 than in the Milky Way or the V-I 
color is systematically 0.09 mag redder. The reason for this difference is unclear. We found for the 
three nights (the subscripts denote the nights): 



®http://www.eso.org/observing/dfo/quality/ FORS/qc/zeropoints/zeropoints.html 
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Vi 
Ii 

V2 
I2 

V3 
I3 



= Vinst + 28.057 + (0.74) - 0.16Xv 
iinst + 27.284 + (0.74) - O.I6X1 - 0.07(Vi 



Ri) 



= Vinst + 28.048 + (0.74) - O.lGXy 
iinst + 27.263 + (0.74) - O.I6X1 - 0.07(V2 



R2) 



= Vinst + 28.038 + (0.74) - 0.16Xv 
iinst + 27.261 + (0.74) - O.I6X1 - 0.07(V3 



R3) 



We note the term due to the conversion from ADUs to electrons (0.74) separately to facilitate 
the comparison with the zero points given in the ESO Web page. 

For the reddening towards NGC 1399 we adopted, according to Schlegel et al. ((1998)), Eb-v = 
0.013. Using Ec-Ti = 1.97 Er-v (Harris & Canterna (1977)) we have to correct C-Tl by 0.026 mag. 
For the comparison field the foreground reddening is negligible: Eb-v = 0.002. 

The color magnitude diagrams for the point sources in the MOSAIC images are plotted in 
Fig. 3. The combined color-magnitude diagram for the three VLT fields is shown in Fig. 5. 



The globular cluster luminosity function (GCLF) will be derived from the VLT images only, so 
the photometric incompleteness as a function of apparent magnitude of the MOSAIC images is not 
of much interest to our current study. However, it is interesting to know the global incompleteness, 
and how many objects, being resolved in the VLT images, were falsely classified as point sources 
in the MOSAIC images. We found that within our magnitude range, (41 it 3)% of the cluster 
candidates visible on the VLT frames were detected on the MOSAIC images. 0.5 objccts/arcmin^ 
which are resolved in the VLT images have been classified as point sources in the MOSAIC images 
(objects brighter than the employed faint magnitude limit). This is a reassuringly small number. 

Dithering to fill the gaps between the eight individual CCD chips was done only for the C 
images. Therefore one has to consider the spatial incompleteness due to the gaps between the 
MOSACI CCDs in the R images, saturated galaxy centers and some bright stars west of NGC 1399 
which account for a spatial incompleteness of Ri 4% in the MOSAIC data. 

Experiments with artificial stars were performed in order to check on whether a magnitude- 
dependent difference between input and output magnitude is observed among the point sources. 
This was not the case, and thus the size of the median filter we used to model and subtract the 
galaxy light, is justified. Additionally, we found that for point sources more distant from NGC 1399 
than 1.6', no radially varying completeness has to be corrected for. 

The completeness of the VLT fields has been determined with extensive artificial star tests. We 



2.4. Photometric Completeness 
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used addstar of the DAOPhot package under IRAF to add 5 times 2000 stars to each VLT image. 

The stellar PSF was the same as determined by the PSF fitting routine in the same package. The 
color and magnitude range encompassed the range of observed globular clusters. The resulting 
completeness as a function of V magnitude is shown in Fig. 6. 

3. The Color Magnitude Diagram 

We will now discuss the CMD presented in Fig. 3 and Fig. 5. The globular clusters show 
up nicely in Fig. 3 within the color range 1.1 < C — Tl < 2.3. In its vast majority the cluster 
population belongs to NGC 1399, but the GCSs of NGC 1404 and NGC 1387 are also included. 
The combination of a large area, depth, and broad color baseline render this CMD unique among 
the available photometric data sets of cluster systems. The cluster population starts at about 
Tl=20. Well visible is the blue "peak" in the color distribution, while the red clusters rather build 
a "shoulder" and do not distinctly peak at a certain color (see section 5 for a detailed discussion 
of the color distribution). The comparison with the background field (Fig. 4) shows how well it 
reproduces the shoulder starting at C-Tl = and the red peak (stars) at C-Tl = 3. 

The comparison between the VLT and the MOSAIC data also shows why the contamination 
of the cluster system by background galaxies is much lower for the filter combination C,R than for 
V,I, even if the seeing is considerably worse. The reason can be seen in Fig. 7. The main sources 
that contaminate the cluster sample in V & I are unresolved faint blue galaxies that have an UV 
excess and are thus bluer than cluster candidates in C-R while having the same V-I color. We 
confirm the galaxy nature of these (C-Tl)-blue objects since they have a homogeneous distribution 
and are not concentrated towards NGC 1399 like the cluster candidates. 

4. The Globular Cluster Luminosity Function 

To derive the globular cluster luminosity function (GCLF) of the VLT data we selected clus- 
ters within a radial distance of 4.3' from the center of NGC 1399. As a background field we used 
the F0RS2 image centered 13' away from NGC1399. We are aware of the fact that this "back- 
ground field" still contains a considerable number of clusters apart from background objects, which 
deteriorates the statistics, but should not alter the location of the turn-over magnitude (TOM). 

The luminosity function is also plotted as kernel filtered density distributions. We followed here 
the approach described by Merritt & Tremblay ((1994)) and used an adaptive Epanechikov kernel. 
The scaling in Fig. 8 is number of clusters per 0.2 magnitude. In the following, all distributions using 
this nonparametric estimation are based on the Epanechikov kernel. For the fitting we use both, 
filtered data and the histograms to ensure that the results do not depend on the representation. 

We fitted a Gaussian to the distributions for different limiting magnitudes (always brighter 
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than the 50% completeness) and found an average a = 1.30 it 0.12 and niTO = 24.0 it 0.1. However, 
a correlation between a and TOM exist these fits which is caused by the fact that our observations 
do not extend far enough beyond the TOM. It is therefore advisable to keep the width fixed and to 
fit only the TOM. With a a = 1.2, 1.3, 1.4, which are used for elliptical galaxies (Ashman & Zepf 
(1998)), we find a TOM of 23.98 ± 0.11, 24.02 ± 0.09, 24.05 ± 0.12. We choose the same final TOM 
of mTO = 24.0 ± 0.1. 

With the knowledge of the TOM we calculate that the fraction of all clusters which are found 
in the area covered by the VLT data is 85% it 4%. 

The background-subtracted luminosity functions of the MOSAIC cluster candidates from 2' to 
16' (excluding 2' around NGC 1404) are shown in Fig.9. The luminosity distributions are consistent 
with the assumed distance modulus of m— M = 31.4 and a width of the Gaussian luminosity function 
of 1.2. 

This TOM together with the known completeness of the VLT data can be used to determine 

the completeness of the MOSAIC fields. Since 41% it 3% of the VLT candidates were found on the 
MOSAIC images (see above) the absolute completeness is 36% ± 5% for the MOSAIC images. 

An alternative approach to find the completeness is to use the luminosity function of the 
MOSAIC clusters directly: using a TOM of Tl= 23.3 and a width of o" = 1.2 (Forte ct al. (2001)) 
we find a completeness of 41 ±2%. We adopt a final completeness of 34% it 5%, taking into account 
a spatial incompleteness due to the undithered gaps of 4%. 

The GCLF is frequently used as a distance indicator for early-type galaxies due to its universal 
peak luminosity (e.g. Harris (1991), Harris (2001), Ashman et al. (1995), Kundu & Whitmore 
(2001a), Larsen et al. (2001)). However, due to the uncertainty of the VLT calibration we do not 
use the obtained TOM for more than the above presented completeness calculations for which the 
absolute calibration is not critical. 

In the following we use a distance to NGC 1399 of m - M = 31.4 ± 0.2 consistent with SBF 
distances (Tonry et al. (2001), Liu et al. (2002)), however only marginally consistent with the 
result of Grillmair et al. ((1999)), from the GCLF technique. 
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5. 



Color Distribution of the Clusters 



In the CMD of all point sources around NGC 1399 (Fig. 3) the GC candidates stand out 
strikingly. A blue peak at C-Tl = 1.3 followed by a broad distribution towards redder colors is 
prominently visible. In the following, we define GC candidates by the magnitude interval 20 < 
Tl < 23 and the color interval 0.8 < C - TK 2.3. 

In Fig. 10 we plot the background subtracted color distribution for the cluster candidates for 
the entire sample and in three different ranges in radial distance. It is apparent that the color 
distribution is radially dependent: at larger radii the relative contribution of the red clusters is 
much smaller than at lower radii. This already shows that the blue clusters exhibit a shallower 
density profile than the red clusters, as has already been found by several authors. However, it is also 
apparent that this might not hold for larger radii: the two outer distributions (the two lower panels 
in Fig. 10) look quite similar and are statistically not distinguishable. The distributions plotted in 
Fig. 10 are tabulated in Table 2. The two peaks in the innermost sample are at C — Tl = 1.32 ±0.05 
and 1.79 it 0.03, respectively. 

Since so far only a few GCS have been studied in the Washington system, our possibilities for 
comparison with other work are quite limited. 

As mentioned above, NGC 1399 has already been studied by Ostrov et al. ((1998)) who found 
a bimodal distribution with peaks at C-Tlf« 1.4 and C-Tlw 1.8, agreeing with our values for the 
inner region. This is also shown in the upper right panel of Fig. 10. 



The area we cover is the largest ever observed for the NGC 1399 GCS. However, the innermost 
region is unaccessible due to the brightness of the galaxy and the saturated parts of the images. 
Therefore we combined our data with HST observations by Forbes et al. ((1998)) in order to 
include also the inner part of NGC 1399 in our study. We extrapolated all three samples (HST, 
VLT, Mosaic) to the entire luminosity function with the parameters given in section 3. In Fig. 12 the 
resulting surface densities of GCs are plotted for the three samples. The HST results are in excellent 
agreement with our data at overlapping radii, demonstrating the validity of our completeness 
calculations. A fit which satisfactorily reproduces the surface densities of the total number of 
clusters after completeness correction is: 



In the following, wc focus on the MOSAIC sample. From the color histograms in Fig. 10 one can 
immediately sec that the red clusters are more centrally concentrated than the blue ones. Guided 
by Fig. 10, we define the color limit between red and blue clusters to be C-Tl = 1.55, which is the 



6. 



Radial Distribution of the Clusters 



n(r) = 200 • (l + {r/l' f) [arcmin~2] 



(1) 
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minimum seen in the clearly bimodal color distribution seen in the upper right panel of Fig. 10. 

We can now quantify the radial surface density distributions displayed in Fig. 13, which show the 
density profiles for both populations. Because the inner region, where a core is visible, is excluded, 
we can describe the density profiles by pure power-laws. 

Out to a galactoccntric radius of about 7', the blue clusters show a distinctly shallower slope 
than the red clusters. This has already been found in earlier studies of NGC 1399 (e.g. Forbes et 
al. (1998), Ostrov et al. (1998)) as well as in other galaxies, NGC 4636 being the most striking 
example (Kissler et al. (1994)). 

The region beyond 7' has as yet been unexplored and indeed we see an interesting change. 
The distribution of the red population is well represented by a single power-law over the full radius 
interval between 2' and 20', but the blue population shows a change from a shallow profile in the 
inner region to a steeper profile in the outer region. Here, two power-laws are more appropriate. 
For the blue clusters, the profiles are proportional to /.-O.SiO.i? ^^^^ j,-i.7i±0.2i ^^^j. ^ 5' < j- < 7' 
and 7' < r < 21', respectively). For the red population, we find the profile ex 7--i-64±o.io_ Thus 
we cannot see any significant difference in the slope between red and blue cluster distributions 
beyond 7' of radial distance. This can also be quantified with a K-S test for which we selected 
only clusters brighter than R=22.5 to minimize the effect of the background. This test results in 
a 48% probability that the red and blue clusters at large radii have a different radial distribution. 
This probability should be considered as upper limit, since we did not correct for the background 
contamination which affects the red and blue color range differently. This indicates again that the 
distribution of red and blue clusters is consistent with being the same. 

One might recognize that the slopes fitted to the MOSAIC data at large radii arc shallower 
than the fit to all GCs over the whole radial range presented above. This is also apparent from 
Fig. 12, which suggests that the assumption of a uniform power-law over the whole galactocentric 
distance range might not be adequate. At larger radii - where the investigation is restricted to the 
MOSAIC data - the slope is shallower. 

To provide an independent check, we derived the radial profiles for the red and blue subsam- 
plcs between 1' and 8' also from the VLT data. We only used cluster candidates brighter than 
V=22 mag to avoid the background contamination coming from the faint blue galaxies, which in 
V-I cannot be distinguished from cluster candidates (sec Fig. 7). Moreover, the VLT fields only 
cover approximately 11% of the MOSAIC field and thus the statistics are poorer and the errors 
larger. In spite of this, the resulting density distributions of metal-poor (0.7 < {V — I) < 0.95) 
and metal-rich clusters (1.05 < {V — I) < 1.2) agree very well with the results obtained with the 
Washington observations: The exponents are —1.0 it 0.13 and —1.5 ± 0.18 for the blue and the 
red populations, respectively. However, a difference emerges for the blue population at small radii. 
If we only consider the blue clusters within 5', we obtain an exponent of —0.5 ± 0.09 compared 
to — 1.2 ± 0.2 (based on 4 points) outside this radius. This indicates that the assumption of two 
power-laws describing the metal-poor population is probably too simple. 
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Besides NGC 1399, NGC 4472 - the brightest galaxy in the Virgo cluster - is the only elliptical 
galaxy where the GCS has been studied with wide-field CCD photometry (Rhode &:Zepf (2001)). 
Due to its similar distance (17 Mpc) and similar size the results for this galaxy can readily compared 
with ours. In particular the radial distributions of red and blue clusters found by Rhode &:Zcpf 
((2001)) is very similar to our results for NGC 1399. This includes the mentioned "out-leveling" of 
the number ratio of red versus blue clusters at large radii, i.e. the same radial decline of red and 
blue clusters is indicated. 

6.1. The Brightest Clusters 

In Fig. 11 we show the color distribution for three different brightness bins. It is striking that 
for the brightest clusters, shown in the uppermost panel, the color distribution differs largely from 
the fainter samples seen (this is already apparent in the CMD): compared to the fainter sample 
mor intermediate colored (C — Tl ~ 1.55) are present. A Kolmogorov-Smirnov test (K-S test) 
resulted in only a 5% probability that the bright sample has the same color distribution as the 
fainter one (R=21-22). Ostrov et al. ((1998)) suggested that the brightest clusters might have a 
unimodal distribution, however their lower number of objects did not permit a clear statement. In 
the faintest magnitude bin, the red peak is shifted towards the blue with respect to the intermediate 
bin. This, however, may be easily explained by the inclined completeness line that docs affect only 
this bin. These bright clusters are as concentrated as the red clusters and do not have the shallow 
distribution of the blue ones: a K-S test results in a 1% probability that the bright ones are not 
distributed than the red ones and a 95% probability that they have a different distribution than 
the blue ones. 



7. The Galaxy Light Profile 

It is now of interest to compare the surface density profile of the GCS with that of the galaxy 
light and to search for any color gradient, which would indicate a correspondence between the GCS 
and the galaxy population itself. 

Previous photographic work covering a similarly large field has been carried out by Schombert 
((1986); V-band) and Caon et al.((1994); B-band). The resulting profiles are plotted in Fig. 14. 
Beyond a radius of 5' they begin to diverge in the sense that the V profile gets brighter than the B 
profile. Besides the fact that this "color gradient" has the very unexpected sense of getting redder, 
it also reaches the unplausibly high value of almost 2 mag at 14'. We conclude that either V or B, 
or both, are not reliable at faint surface brightness values. We therefore tried to derive the light 
profile from our images as well. 

We measured the brightness around all point sources within an annulus between 8" and 11" 
on the point source-subtracted C and R images. We defined this brightness to be the mode of the 
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distribution of pixel intensities after a 5% clipping of the brightest peaks. Without the a clipping 
the scatter would have been significantly larger. 

Care has to be taken due to the use of eight different chips and corresponding differences in 
the sensitivities that arc left after flatfielding. To account for these differences we shifted the zero 
points of individual chips to a common mean value by adjusting the sky brightness in adjacent 
fields. The resulting correction factors are of the order of 0.3% of the background level in the R 
image and in the order of 1% in the C images. Since the R images are deeper we used them to trace 
the galaxy light to large radii, whereas the C data has been used to search for a color gradient at 
smaller radii. We excluded regions with bright stars, other galaxies and obvious fiatfield structure 
(chip #4, chip #5 and upper part of chip #8). The inner 2' are saturated. 

The most difficult part in determining the galaxy luminosity profile is the correct determination 
of the sky background: a tiny 0.5% deviation in the sky level has already a profound effect on the 
galaxy light distribution as shown in Fig. 14 for the Tl profile. In the upper panel the luminosity 
profile is given for different sky values; these sky levels are indicated in a linear intensity scale in the 
lower panel, where the x-axis is enlarged. Individual measurements are plotted as dots. If we would 
have reached the true sky background, profile b would be appropriate and our resulting luminosity 
profile would look nearly identical to that of Caon etal. ((1994)). However, the measurements 
suggest that the galaxy light has not yet leveled out. On the other hand profile a can be excluded, 
since the observed background is clearly lower that the one belonging to this profile. The profile 
by Schombert et al. can be reproduced (profile d) when using the lowest background in the lower 
panel of Fig. 14, but such a low background is almost certainly excluded since in this case the 
galaxy would have a radial extension of approximately 500 kpc and hence it would be one of the 
largest galaxies known. As will be discussed later on, profile c, based on our prefered baclground 
also resembles well that of the GCs and thus wc adopt it as our preferred profile. The difference 
with the B-profile of Caon et al. would still suggest a color which becomes redder with increasing 
distance, contrary to what we find (see below). Only profile a delivers a color gradient in the 
required sense, but this is not supported by our data. 

From Fig. 14 it is also apparent that the profile becomes rather diffuse for radii smaller than 
approximately 6'. The reason for this is the non-negligible ellipticity of NGC 1399. The lower 
profile envelope belongs to the minor axis while the upper one belongs to the major axis. We 
determined an average decline only, but this poses no problem for the comparison with the cluster 

distribution, which has also been determined assuming a circular geometry. Moreover, the slope of 
a power-law is not altered by averaging over the ellipticity. We nevertheless used this information 
to derive the radial dependence of the ellipticity which is shown in the next section. 

The C-profilc reaches the background level already at 13' due to the shallower C-data and we 
determine the background by averaging the "profile" at larger radii. 

Between 2.5' and 9' (where the effects of the background uncertainty are still tolerable) we 
fitted power-laws to the profile and obtained exponents of — 1.85 ± 0.03 for the Tl profile and of 
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—1.21 lb 0.02 for the C profile, thus finding a clearly shallower C-profile. This holds for any sky 

background which deviates no more than 2% from our measured background, to which we assign 
an error of less than 1%. The color gradient in C-Tl is shown in the lower panel of Fig. 20 in 
Sec. 7. 2. The color gradient is also plotted as a smoothed curve, together with error limits that are 
estimated from the background uncertainties, in Fig. 21. 

We see that the color gradient might level out to a constant value of C-Tl around 1.2 at 
6', which is approximately the same radius beyond which the blue and red GCs start to show 
indistinguishable surface density profiles. This might be evidence that the behavior of the GCS is 
indeed reflected in the field star population. However, the mean color of the clusters at this radius 
is (linearly averaged) C-T1=1.47, and hence considerably redder than the galaxy color. 

7.1. Covering the whole Radial Range 

We are unable to determine the total magnitude of NGC 1399 because of the missing inner 2'. 
However, we can proceed by taking this data from the literature. In Fig. 15 we plot the available 
luminosity profiles which were shifted onto our Tl-magnitude scale using our Tl-profile. We shifted 
Caon et al.'s and Schombert et al.'s data in the range between 1' and 1.5'. In the innermost region 
a V-profile based on HST observations has been published by Lauer et al. (1995), which we used 
out to 12" . In the range between 12" and 1.5', we used the B profile of Caon et al. (1994). For 
larger radii, we used our profile. Color gradients in B-R and V-R in the inner region are very small 
and should not introduce significant errors (e.g. Idiart et al. (2002)). 

From this data we can construct a luminosity profile for NCC 1399 that covers 3 orders of 
magnitude in radius. In Tab. 10 we list this profile, together with the integrated Tl magnitude at 
selected radii. 

A good fit to the surface brightness profile of NGC 1399 is: 

Tl{r) = 2.125 log(l + (r/O.Obb'f) + 15.75. (2) 

The deviation from the fit is shown in Fig. 16, where we plot the difference between observed 
and fitted magnitude. Similarly to the globular cluster density profile, the fit to the light profile 
using the whole radial range also results in a profile being too steep at large radii (see Fig. 16). For 
these large radii, an exponent of 1.83 as derived above is more adequate. Integrating this profile, 
one obtains the total R-luminosity of NGC 1399: Mr = —23.33 (assuming a distance modulus of 
31.4 and an absolute solar luminosity of Mr^ = 4.28.). 

If we adopt spherical geometry, one can derive the deprojcctcd profile. Because we have the 
profile as an analytic expression, an easier way than to apply the dcprojection formula is to adopt 
a power-law of the above form, use the projection formula and vary the parameters until the best 
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fit is found. 

Using this procedure, the best fit we obtain is: 

L{T1)[Lq/pc^] = 101 • (1 + r/221pc)-2-85. (3) 

We also determined the ellipticity of NGC 1399. However, we did not try to fit a position angle 
to the rather noisy data, but estimated the position angle rather by eye inspection of the image. 
We used a fixed position angle of 90°. The position angle variation from 110° to 85° seen by Caon 
et al. ((1994)) is not critical since it occurs between 30" and 1.35', while our profile starts at 2'. 
The resulting radial dependence of the ellipticity is shown in Fig. 17 together with Caon et al.'s 
values. The two measurements are in good agreement. The ellipticity of NGC 1399 is 0.1 between 
2' and 4', then it rises to about 0.2 at 5'. The further decline at larger radii is rather uncertain, 
since at these large radii our azimuthal coverage is not homogeneous due to the exclusion of areas 
with bright stars, other galaxies and obviously remaining flatfield structure. 

In Fig. 18 we compare the integrated light profile with aperture measurements from the litera- 
ture. This figure shows that our determination agrees well with the literature values, and that our 
data do not suggest any systematic zero point difference. 

8. The Specific Frequency 

The specific frequency is the number of globular clusters per unit luminosity scaled to a galaxy 
with an absolute luminosity of My = -15 mag {Sy = N ■ i{)OA{Mv-ib)y 

In Fig. 12, the galaxy light profile is compared with the cluster density profile over nearly 
the whole observed radial range. Determining the specific frequency means integrating over both 
curves. In the inner part the galaxy light is more concentrated than the globular clusters, which 
reflects the fact that the core radius of the GCS is about 18 times larger than that of the field 
population. This also means that the specific frequency is not radially constant. We shall proceed 
in deriving first a global value, and then study the radial dependence of the specific frequency for 
both red and blue clusters. 



8.1. The Global Specific Frequency 

To derive the total number of clusters we have to account for the incompleteness with respect 
to the globular cluster luminosity function. In Sect. 2.4 we determined this to be (35 ± 5)%. 

Within the radial range from 2' to 17' (excluding a radius of 2' around NGC 1404), we find 
2600 ±50 cluster candidates. Within the central 2' Forbes et al. ((1998)) found 700 clusters. From 
our background field we know that a cluster candidate sample brighter than Tl = 23.5 has a con- 
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tamination of 0.48 it 0.02 objects/arcmin^. Accounting for this contamination and incompleteness, 
we find a total number of 6100 zt 770 clusters within 83kpc (15'). 

Alternatively, we integrated the density profile derived in Sec. 7 and obtained a number of 
6800 ±950 clusters within 83kpc. Within 10', this results in 5793 clusters. Averaging the results of 
both determinations we find the total number of clusters within 15' to be 6450 it 700. Our obtained 
number agrees well with what has been found for the cluster content within 10': Ostrov et al. 
((1998)), Forbes et al. ((1998)), and Kissler-Patig et al. ((1997)) all quote ^ 5700 clusters within 
10'. In the inner part we have an excellent agreement with older studies, too: within 1.5' and 2.5' 
we determine a total number of 734 ± 114 clusters, while Ostrov et al. found 746 it 80, and within 
2.5' and 4' we have 921 ± 140 whereas Ostrov et al. give 967 ± 100. 

With this total number of clusters, we find a global specific frequency of Sr = 3.2 it 0.7 using a 
total galaxy magnitude of Tl= 8.08 (Table 3) and m-M=31.4 it 0.2. Since we cannot measure the 
galaxy light in the V band, we have to convert Tl into V magnitudes to be able to compare our 
value with the literature. Using a color of V-R= 0.55 (Mackie (1990)) we find a specific frequency 
of Sv = 5.1 ± 1.2 within 15'. 

Forbes et al. ((1998)) quoted a specific frequency of 11.5 it 1. They used an absolute mag- 
nitude of My = —21.74 (taken from Faber et al. (1989)), while we would use Tl=8.25, which 
translates to My = —22.59 (see above) for the total luminosity (within 10'). Part of the difference 
comes from the different distance modulus adopted by Forbes et al. (m-M=31.2). However, the 
major discrepancy is in the apparent magnitude: Faber et al. give the Bt, magnitude which is 
defined as the magnitude within the 25 mag/arcsec^ isophote. Using a B-R=1.8 value, we find that 
25 mag/arcscc^ corresponds to a radius of 3.8'. Within this radius we measure an approximate B 
luminosity (again assuming no B-R color gradient) of 10.4 mag, which is in reasonable agreement, 
considering all uncertainties to the Faber et al. value of B=10.55. The same argument can also be 
applied to explain the high specific frequencies determined by Kissler-Patig et al. ((1997)), Wagner 
et al. ((1991)), and Bridges et al. ((1991)). This shows that one of the most important caveats in 
deriving the specific frequency is that the clusters and the galaxy light has to be measured within 
the same area, a fact that has already been pointed out by Ostrov et al. ((1998)) for the case of 
NGC 1399. 

8.2. The Local Specific Frequency 

Because of the existence of subpopulations in the GCS, it is useful to define a local specific 
frequency where the number of GCs is normalized to the brightness of the area in which they are 
counted. 

The local specific frequency of the red, the blue and the total cluster sample for different 
radial bins is given in Tab.4 and displayed in Fig. 19. Since the galaxy light profile in the R 
band matches well the radial distribution of the red cluster population within 11', its local specific 
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frequency remains nearly constant around Sr = 3. Outside of 11' small background uncertainties 
make a definite statement on the surface brightness profile difficult. However, using a reasonable 
background it is easy to find a galaxy profile that follows the red clusters also at these large 
radii. The blue cluster system, on the other hand, shows within 11' a slightly shallower profile 
and therefore its specific frequency with respect to the R band luminosity increases outwards 
proportional to 7-0-8±o.2 jjowever, the behavior at larger radii is uncertain as well. 

9. Discussion 

9.1. The Cluster Color Distribution 

In general, there is little doubt that the color of a GC indicates its metallicity in an old elliptical 
galaxy like NGC 1399 (however, we cannot exclude that younger clusters are mixed in; Forbes et 
al. (2001b)). The GC color distribution has been the starting point for many investigations 
regarding the substructure of GCSs of elliptical galaxies. It has been found by many authors (e.g. 
Kundu & Whitmore (2001a), Forbes & Forte (2001a), Larsen et al. (2001), Gcbhardt & Kissler- 
Patig (1999), see also the corresponding section in the book of Ashman and Zepf (1998)) that 
GCSs often exhibit two peaks in the color distribution, labeled as "bimodality" . If bimodality 
is not obvious then frequently at least the existence of two peaks is statistically better supported 
than an unimodal color distribution. Usually, a linear transformation between color and metallicity 
without accounting to an intrinsic scatter around the mean color metallicity is applied (e.g. Harris 
et al. (2000), Kissler-Patig et al. (1997), Harris et al. (1991)). This means that a bimodal color 
distribution translates into a bimodal metallicity distribution which is assigned to two different 
populations. Because bimodal color distributions of clusters in ellipticals have been predicted by 
Ashman &; Zepf (1992) in the context of a simple merger scenario, color bimodality has been 
interpreted supporting the formation of GCSs of elliptical host galaxies in merger events. However, 
despite the fact that a merger can produce a bimodal color distribution it is not clear whether 
a bimodal distribution necessarily implies a merger scenario, e.g. even in a monolithic collapse 
scenario a bimodal metallicity distribution (of stars) can be obtained: Samland et al. ((1997)) 
simulated the formation of a disk galaxy starting with a coherent gas cloud. They found that 
disk, halo and bulge stars have different metallicity distributions (their Fig. 7). If one plots their 
Fig.7 in logarithmic metallicity bins a bi- or even multimodal metallicity distribution becomes 
apparent. They note "that a bimodal star formation or a pre-enrichment of the protogalactic cloud 
is not required by our model. The metallicity distributions for the different galactic regions result 
pTirely from the local star formation history, the large-scale gas flows, and the phase transitions by 
condensation and evaporation." 

We also would like to call attention to the fact that the color-metallicity relation is non-linear, 
flattening out in the metal-poor regime (Harris & Harris (2002)). In conjunction with an intrinsic 
spread in the color-metallicity relation due to second parameters and photometric erros, this leads 
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to the circumstance that, starting from the bluest colors equidistant color intervals are projected 
onto progressively larger metallicity intervals, which decrease once the linear regime is reached. 
What this effect means for the interpretation in particular of the blue peak, is certainly worth to 
investigate, but beyond the scope of our paper. 

The question whether the color distribution of a GCS depends on the host galaxy's properties 
is still unclear. Regarding the Washington system only few GCSs have been studied: Geisler et 
al. ((1996)) investigated the GCS of NGC 4472 and obtained C - Tl ?a 1.3 and C - Tl « 1.8 for 
the colors of the blue and read peak, respectively. The blue peak of NGC 1427, a low luminosity 
elliptical galaxy, also agrees with the blue peak of the NGC 1399 clusters (C — Tl « 1.4, Forte et 
al. (2001)). M87 also has a GCS with a bimodal color distribution that peaks at C-Tl?a 1.25, 
while the position of the red peak remains uncertain in the color histogram shown by Cote et al. 
(2001), which is based on still unpublished photometry. Two more GCSs for which Washington 
photometry have been obtained are those of NGC 3311 (Seeker et al. (1995)) and NGC 3923 (Zepf 
et al. (1995)). Both galaxies (which have been observed during the same run) show abnormally 
red color distributions. However, Brodie et al. (2000) performed V-I photometry with HST in 
NGC 3311 and found a normal bimodal distribution with peaks at the expected colors. Since the 
Washington run was not photometric, one might conjecture that both NGC 3311 and NGC 3923 
were subject to calibration uncertainties. Together with our results on NGC 1399 we conclude that 
at least the blue peaks have very similar colors in those GCSs of elliptical galaxies investigated so 
far with Washington photometry. 



9.2. Interpreting the Radial Cluster Distribution 

In Paper II, we will show that the blue and red clusters within 7' are also kinematically dis- 
tinct: The red clusters have a velocity dispersion of 274 it 15km/s, while the blue clusters have 
310 ± lOkm/s. This larger dispersion can well be accounted for in the framework of a spheri- 
cally dynamical model due to the different surface density slopes. In other words, a higher/lower 
dispersion is dynamically equivalent to a shallower/steeper surface density profile in an isotropic 
equilibrium situation. 

A distinction of two subpopulations on the basis of the color histogram alone bears a high 
degree of uncertainty. However, in light of the correlation between metallicity and density pro- 
file/kinematics, we find it justified to talk about two populations out to a radius of about 7'. 
However, that does not necessarily imply the existence of two distinct popTilations, in the sense of 
different formation epochs or formation mechanisms. It may well be that there is a continuum of 
properties, over which our color sampling is simply taking an average. Going to larger galactocentric 
radii, the distinction between subpopulations becomes much less pronounced. The density profiles 
of blue and red clusters are more or less the same and the color distribution indicates a broad range 
of metallicities, even if two Gaussians are still a better representation that one Gaussian, but as we 
said above, we do not assign too much physical significance to that. 
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Despite the large field used we have not yet reached the background. The total extension of the 
GCS of NGC 1399 is hence not yet known. However, for GCs with a very large radii the question 
remains whether such clusters should be considered as belonging to NGC 1399 or as an " intergalactic 
GC population" (White (1987), West et al. (1995)). Bassino et al. ((2002)) used Washington 
observations around dwarf elliptical galaxies in the Fornax clusters to search for GCs at large radii. 
They found a GC candidate surface density of 0.25 it 0.08/arcmin^ and 0.12 it 0.06/arcmin^ at a 
radial distance of 40' and 2° from NGC 1399, respectively. For comparison, using their limiting 
magnitude of R=22mag, we find a surface density of 0.3 it 0.1/arcmin^ at a radial distance of 16'. 
Extrapolating our derived radial profile we would expect a surface density of ~ 0.14/arcmin^ at 40' 
and of ~ 0.02/arcmin^ at 2°. In particular at the larger radius we would expect a smaller number 
of GCs. This could be an indication that a population of "intergalactic clusters" exists, as proposed 
by Bassino et al. ((2002)). 



9.3. Comparison of the NGC 1399 Light Profile with that of a "Normal" Elliptical 

The luminosity profiles of many "normal" ellipticals follow a r^/^-de Vaucouleurs law. Here 

we have the opportunity to differentially compare NGC 1399 with such a "normal" galaxy, NGC 
1404. This difference is illustrated in Fig. 20, where the profiles of NGC 1399 and of NGC 1404 
arc compared. Wc derived the luminosity profile of NGC 1404 in C and Tl in the same way as 
for NGC 1399. Care has been taken to take account of the extended halo of NGC 1399, which 
envelopes NGC 1404 fully in projection. We subtracted its contribution, i.e. the final luminosity 
profile, from each pixel. For the sky background the same value as for NGC 1399 has been used. 

It can be seen that the light profile of NGC 1404, in contrast to NGC 1399, does not exhibit a 
uniform power law, but has an exponent of about -2 for radii smaller than 1.8' and -2.9 for larger 
radii, i.e. almost -4 in deprojection. 

The deviation from a r^/^-de Vaucouleurs law is not unusual for galaxies brighter than My = 
—21. The light profile shapes of brighter galaxies become progressively more concave when plotted 
against r^^^ (Kormendy & Djorgovski (1989)). These authors define a cD galaxy by a light profile, 
which exhibits an inflection point independent from the way it is plotted. In this sense, NGC 1399 
would not be a cD galaxy. 

Whatever the exact definition of a cD galaxy may be, it is useful for further insight to think 
about the physical reason of a de Vaucouleurs profile. Many "normal" elliptical galaxies can be well 
described by a Jaffe (1987) profile, which in projection resembles the de Vaucouleurs law, and which 
in space in the outer regions is characterized by a transition from a to a power law. Jaffe 
((1987)) and White ((1987)) have shown that a power law emerges if the energy distribution 
function has a sharp break at the escape energy, given a nearly Keplerian potential. 

NGC 1399 has a massive dark halo resembling the potential of an isothermal sphere out 
to a distance of at least 40 kpc, shown by the dynamics of the globular clusters (Paper ll). The 
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temperature of the X-ray gas outside this radius remains isothermal (Jones et al. (1997)). Therefore, 
even if the exact shape of the potential is not known, we can at least conclude that the Keplerian 
regime is not reached out to the very faint surface brightnesses where our photometry becomes 
unreliable. Thus it is very tempting to relate the persistence of the profile to the dark halo of 
NGC 1399. 

We therefore understand the "halo" of NGC 1399, which morphologically appears as a deviation 
from a de Vaucouleurs profile, not as an additional component but as the natural continuation of 
the galaxy profile in the isothermal potential of the dark matter halo. Thus the term "cD halo" 
might be misleading since it implies an additional stellar component. 

Finally, we do not detect a color gradient for NGC 1404. This shows that the population 
composition of NGC 1404 is radially more homogeneous than that of NGC 1399. Additionally, this 
also demonstrates that the observed color gradient in NGC 1399 is not an artefact of an incorrect 
background subtraction. 

9.4. Disentangling the Stellar Populations 

Superior to the approach to relate the two observed red and blue cluster subpopulations to the 
integral galaxy light would be the comparison between clusters and the respective field populations 
that formed during the same star formation event. This is not possible in a direct way. However, 
we can follow a more indirect approach in which we model the galaxy C and R light profile using 
some simplifying assumptions. We assume the field population to be a superposition of only two 
populations, one corresponding to the blue clusters and one to the red clusters. These two field 
populations shall have the mean color of their corresponding cluster populations. The relation 
between the number of clusters and the luminosity of their corresponding field population will be 
called intrinsic specific frequency (ISF) and is measured with respect to the Tl filter. 

It is not possible to fit both, the color and the luminosity profile, if both cluster populations 
would have the same constant ISF because of their different surface density profiles. In addition, to 
reproduce the color gradient requires to allow for a change of the properties of the red populations by 
either a radial color or radial local specific frequency dependence. We limit ourselves to consider a 
changing ISF for the red population, which might be the strongest assumption we make. Otherwise 
the problem would not be well constrained and many solutions could be found. We used the linear 
approximation ISF oc r. A shallower variation, e.g. ISF oc log(r) results in a too shallow color 
gradient. With a much steeper dependence, e.g. ISF oc r^ the model luminosity profile shows a 
considerable dip at the point where the light of the blue population begins to dominate the total 
light, which is also incompatible with the observations. A model that fits the data reasonably well 
is depicted in Fig. 21. A characteristic of any satisfactory model is that the ISF of the blue cluster 
population is much higher than that of the red one at smaller radii. 

It is also possible to use the color and luminosity profile together with the density profiles of red 
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and blue clusters to derive the ISF for both clusters samples directly. However, the disadvantage 

of this approach is that it is much stronger influenced by noise. Nevertheless, it may serve as a 
consistency test and we include the results of this calculation in Fig. 21. This illustrates that the 
two methods agree well considering the uncertainties. 

9.5. From Light Distribution to Mass Distribution 

We note that the color gradient in NGC 1399 indicates a radial change in the stellar M/L 
profile. If the outer regions are more metal-poor the M/L value decreases towards larger radii. 
If one wants to use the luminosity profile to infer the mass profile, one has to account for this 
effect. For an estimation, we assume two populations with the same age of 12 Gyr of different 
mean metallicities: a metallicity of -0.4 dex (representing the inner part in accordance with the 
galaxy's color) and -1.4 dex (representing the outermost parts). Worthey's interpolation engine 
then gives for the respective M/L-values in R (Miller-Scalo IMF): 1.45 (-0.4 dex) and 1.23 (-1.4 
dex). This means that the associated mass profile becomes somewhat steeper than what would 
be derived from the luminosity profile using a constant M/L. Thus an increase of the power law 
exponent from —1.83 to —2.2 is expected. 

10. Summary and conclusions 

We investigated the GCS of NGC 1399 in the bands Washington C and Kron R with the CTIO 
4-m MOSAIC camera. The imaged area is 36' x 36', thus being among the largest concerning GCS 
investigations. 

The main findings of this study are: 

• The relatively modest specific frequency of 5.1 ± 1.2 of the GCS of NGC 1399, already noted 
by Ostrov et al. ((1998)) is confirmed. The local specific frequency increases with radius and 
is higher for the blue than for the red clusters. 

• The color distribution shows a very pronounced bimodal appearance. However, we argue 
that the blue peak might artificially be caused by the non-linear color metallicity relation in 
conjunction with the color scatter at a given metallicity. 

• The surface density profile is shallower for the blue clusters than for the red clusters within 
8 arcmin. Outside a radius of 8 ', both surface density profiles are not distinguishable. 



The galaxy light shows a colour gradient resembling the change of the surface density profiles 
of blue and red clusters. The galaxy light can be modeled with the simple assumption of 
only two populations, each traced by a cluster population of the same mean color. A model 
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consistent with the data has a constant ISF for the blue clusters and a radially increasing ISF 
for the red clusters. 

• The brightest clusters do not show the pronounced bimodal color distribution of the fainter 
ones and have an at least as concentrated radial distribution as the red GCs. 

NGC 1399 was one of the "classic" high specific frequency galaxies. It would be very interesting 
to perform new studies of the similar cases M87 and NGC 3311 in order to see whether their high 
specific frequency still withstand closer scrutiny. If not, central galaxies in galaxy clusters are 
perhaps not the places of extraordinarily high efficiency of globular cluster formation, for which 
they have been taken. 

The GC color distribution, in particular bimodality alone is not very restrictive for galaxy 
formation scenarios. Various models have been published that explain the observed bimodal- 
ity/metallicity distribution, e.g. Ashman & Zcpf ((1992)) Forbes et al. ((1997)), Cote ct al. 
((2001)), Beasley et al. ((2002)). It is believed that bimodal color distributions can be used as 
an argument against monolithic collapse scenarios, for example Rhode k, Zepf ((2001)). However, 
doubts remain: Samland et al. ((1997)) showed that also a monolithic collapse scenario result in 
populations that have different metallicity distributions. 

Since the color distribution alone is not suitable to constrain formation scenarios it is necessary 
to use further information. 

The close relation of the color gradient of NGC 1399 with the mean color of the GCS might 
be seen in favor of a strongly dissipative period during which the globular clusters and the stellar 

body of the galaxy were formed together. However, a counterargument is the radially increasing 
local specific frequency of the red clusters. During a dissipative collapse one would expect the star 
formation rate to be highest where the density is highest, i.e in the central regions. Accordingly, also 
the efficiency of globular cluster formation should decrease radially outwards, if the star formation 
rate is a determining parameter (Larsen & Richtler (2000)). 

Additionally, one has to take into account the kinematic and dynamical properties of the GCS. 
These are (Paper ii): blue and red clusters show different velocity dispersions in accordance with 
their different surface density profiles. The orbit distribution is nearly isotropic. Only the outer blue 
clusters (if any) show signatures of rotation. The overall mass distribution resembles an isothermal 
sphere. 

Thus some sort of relaxation mechanism must have been active in the early times of NGC 
1399, which is difficult to find in a purely dissipative collapse, but naturally provided by a phase 
of violent relaxation during a merger. The merger components (probably more than only 2) could 
have been gas-rich, unevolvcd metal-poor galaxies where cluster formation was efficient during the 
starburst preceding the merger itself. In addition to the (presumably mainly metal poor) cluster 
populations of the merger components, the bulk of clusters, both metal-rich and metal-poor, could 
have been formed during this period. These clusters have after the merging took place roughly the 
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same spatial distribution. However, after the merger remnant has reached its equihbrium, cluster 
formation continued in the inner regions, albeit with a decreased efficiency, forming the metal rich 
part of the inner red cluster population. 

Based on this idea we can estimate the number of additional red globular clusters. For this we 
assume that we have one population of metal rich clusters that have the same density profile as the 
blue ones (albeit with a different normalization taken from the red clusters at radii larger than 8'). 
Then we require ^ 700 additional red clusters to explain the difference between the two profiles. 
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Table 1. Observational data used in this investigation. 
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Table 2. Color distribution of clusters for the whole sample (Nau), within 1.8' < r < 4.5' (Ni), 

within 4.5' < r < 13.5' (N2) and r > 13.5' (N3). 



C-R 


Nail 


Ni 


N2 


N3 


Nbg 


0.8 


14 ±13 


2±2 


10 ±7 


4±8 


80 ±9 


0.9 


60± 18 


6±3 


31 ±9 


26 ± 11 


129 ± 11 


1 


53± 18 


10 ±4 


41 ± 11 


3±11 


135 ± 12 


1.1 


169 ± 20 


31 ±6 


96± 13 


42 ± 12 


119 ± 11 


1.2 


384 ± 24 


59 ±8 


203 ± 16 


119 ±15 


108 ± 10 


1.3 


458 ± 25 


73 ±9 


253 ± 17 


123 ± 14 


87 ±9 


1.4 


360 ± 22 


72 ±9 


187 ± 15 


88 ± 12 


53 ± 7 


1.5 


223 ± 17 


50 ±7 


109 ± 11 


50 ±9 


42 ± 6 


1.6 


236 ± 17 


55 ±7 


111 ± 11 


55 ± 9 


23 ± 5 


1.7 


186 ± 15 


59 ±8 


83± 10 


37 ± 7 


19 ±4 


1.8 


219 ± 16 


72 ±8 


95 ±10 


41 ±7 


15 ±4 


1.9 


156 ± 13 


50 ±7 


68 ±9 


24 ±6 


10 ±3 


2 


108 ± 12 


37 ±6 


45 ±7 


15 ±5 


13 ±4 


2.1 


39 ±8 


9±3 


18 ±5 


4±3 


9±3 


2.2 


14 ±6 


7± 2 


1 ± 3 


4±4 


12 ± 3 


2.:-! 


5 ± 1 


1 = 1 


:-! ± 2 


1 ± 1 


5 = 2 



-28- 



Table 3. Observed radial number density of red {n(red)) and blue {n(hlue)) clusters. 
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Table 4. The final luminosity profile of NGC 1399 plotted in Fig. 14. The second column gives 

the surface brightness at the radius given in the first column. The third column gives the 
elUpticity and the fourth column the integrated magnitudes at selected radii. The last column 

shows the C-Tl-color. 
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Table 5. Local specific frequency for the whole sample (second column), the red subpopulation 
(third column), and the blue sample (last column). The errors include the error of the luminosity 

profile and Poisson counting statistics. 
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Fig. 1.— 



The MOSAIC field and the three smaller VLT fields are overlayed on a DSS image. 
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Fig. 2. — The residuals of the standard star magnitudes after the apphed caUbration are plotted. 
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Fig. 3.— The CMD based on the MOSAIC data is shown in the left panel for the NGC 1399 field, 
and in the right panel for the background field. The line indicates the limiting magnitude which 
is used throughout the paper. Brighter stars are also plotted with slightly thicker dots than stars 
fainter than the limiting magnitude. 
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Fig. 4. — Comparison of the color histograms of all selected point sources of the NGC 1399 field 
(solid line) and the comparison field (dashed line). 




Fig. 5. — The CMD assembled from the three VLT fields. The bimodality of the color distribution 
is visible, although not as nicely as in the Washington colors (Fig. 3). We shall use the VLT 
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Fig. 6. — Completeness function for the VLT data (Field 1). 
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Fig. 7. — The CMDs of point sources that have been identified on both the VLT and the MOSAIC 
images. We marked objects bluer than C-T1=1 by filled circles. These are predominantly galaxies. 
Note how these objects contaminate the CMD in V,V-I. 
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Fig. 8. — The luminosity function of the cluster candidates within a adial range of 4.3' in the 
F0RS2 images is plotted as solid line. The shaded area designates the magnitude range beyond 
the 50% completeness. The dotted line shows the histogram before the completeness correction. 
The dashed histogram indicates the (arbitrarily scaled) HST GC luminosity function of Grillmair 
et el. ((1999)). 
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Fig. 9. — The luminosity function of the cluster candidates from the MOSAIC data is shown for the 

whole sample (solid line), and the blue and red subsamplcs (dotted and dashed line, respectively). 
The luminosity function given by Ostrov et al. ((1998)) is also shown (arbitrarily scaled). 
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C-Tl C-Tl 



Fig. 10. — The color distribution of globular clusters is plotted for the entire sample and, for three 
samples within different radial intervals which are indicated at the upper right of each panel. The 
dotted histogram shows the distribution before the statistical background subtraction. The smooth 
line is the adaptive kernel filtered distribution (described in Sect. 4). The dashed histogram in the 
upper right panel shows the color distribution by Ostrov et al. ((1998)) within the same area and 
demonstrates the good agreement of the photometry. The counts of Ostrov et al. have not been 
scaled. For an interpretation of these observed distributions see Sect. 5 of the text. 
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C-Tl 

Fig. 11. — The color distribution of globular clusters is plotted for three different luminosity 
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Fig. 12. — This plot shows the radial density distribution of all clusters (extrapolated to the entire 
luminosity function) derived from the MOSAIC data (open circles), the VLT data (filled circles) 
and from HST data (triangles, Forbes et al. (1998)). The solid line is the fit to the data. The 
dashed line shows the luminosity profile of NGC 1399, which is described in Sec. 7. 
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Fig. 13. — The radial density profiles of the observed number of red (filled circles) and blue (open 
circles) globular clusters detected on the MOSAIC image are plotted together with the fitted power- 
laws. The counts are tabulated in Tab. 3. 
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Fig. 14. — The luminosity profile from our R images (solid lines, for four different background 
values), the one published by Schombert et al. ((1986), short dashed line) and the one by Caon et 
al. ((1994), long dashed line). We shifted all profiles to the B surface magnitude given by Caon 
et al. within the radial range of 1' to 1.5'. The lower panel shows the corresponding background 
values for the 4 profiles. 
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Fig. 15. — The Tl-luminosity profile of NGC 1399. We shifted all profiles (even if measured in a 
different passband) to our profile, which is the solid line at large radii. The profile of Caon et al. 
((1994)) is the short dashed curve, the one of Schombert et al. ((1986)) the long dashed curve and 
the HST profile by Lauer et al. ((1995)) is the solid line at small radii. 
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Fig. 16. — Difference between tlie observed and fitted galaxy light profile. Positive values mean 
values brighter than observed. 
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Fig. 17. — The radial dependence of tlie ellipticity of the galaxy Tl light is plotted as solid lines 
in the outer region using our MOSAIC data (using a fixed position angle at 90°) and in the inner 
part using that of Lauer et al. ((1995)). The dashed line shows the ellipticity derived by Caon et 
al. ((1994)). 
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Fig. 18. — Integrated R magnitude dependence on radius for our data and from various literature 
sources. For this purpos the Tl measurements have been transformed into R measurements. 
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Fig. 19. — The radial local specific frequency is plotted for the whole cluster sample and for the 
metal-rich and metal-poor subsamples. 
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Fig. 20. — Upper panel: The luminosity profile in Tl (solid line) and C (dotted line) of NGC 1399 
is be compared to the Tl (dashed line) and C (dash-dotted line) surface luminosity of NGC 1404. 
Lower panel: The radial color is plotted for NGC 1399 (sohd line) and NGC 1404 (dashed line). 
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Fig. 21. — A simple model to fit the galaxy light profile measmed in the C and R band. In the 
upper panel the galaxy light is shown as solid lines accompanied by estimated errors due to 
background variations. The solid circles show the model light distribution, that is a composite of 
the contribution of the red (triangles) and blue (open circles) populations. The two horizontal lines 
indicate the color that has been used as mean color for the two populations. In the middle panel 
the color profile is plotted with a solid line together with the estimated errors and compared to the 
model color distribution (solid circles). In the lowest panel we plotted as open circle and open 
triangle the two ISF that define the model for blue and red population respectively. The crosses and 
the diamonds show the ISF for blue and red populations, respectively, when it is directly calculated. 



